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ABSTRACT 


A  comprehensive  set  of  measurements  has  been  made  in  a  ballistic 
range  which  permits  the  sphere  drag  coefficient  to  be  derived  with  an 
uncertainty  of  approximately  ±2  percent  in  the  flight  regime 
0.  1  <  Ma)<  6.  0  and  2  x  101  <  Re*  <  105  for  Tw/T.  -  1.  0.  Sufficient 
information  is  also  presented  to  predict  the  effect  of  wall  temperature 
on  sphere  drag  coefficient  when  Tw/T00  4  1.0  for  2  <  <  6.  This  in¬ 

vestigation  represents  the  most  comprehensive  experimental  program 
to  date  to  define  sphere  drag  in  the  velocity-altitude  regime  applicable 
to  the  falling  sphere  technique  for  defining  upper  air  density. 
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NOMENCLATURE 


A 

a 

CD 

CDF 

CDPB 

CDPF 

CDFM 

D 

d 

M 

m 

P 

po 

P°inviscid 

R 

Re 

T 

t 

V 

AV 


x 

P 

P 


n 

Frontal  area  of  sphere,  7rd^/4 
Speed  of  sound 
Drag  coefficient 
Friction  drag  coefficient 
Afterbody  pressure  drag  coefficient 
Forebody  pressure  drag  coefficient 
Free-molecular  drag  coefficient 
Drag  force 

Diameter  of  sphere,  (dmax  +  dmjn)/2 

Mach  number,  V/a 

Mass  of  sphere 

Pressure 

Impact  pressure 

Ideal,  inviscid  impact  pressure 

Gas  constant 

Reynolds  number,  pVd/p 

Temperature 

Time 

Velocity 

Velocity  drop  (for  75  ft  of  flight  between  shadowgraph 
stations  1  and  6 

Distance 

Density 

Viscosity 


SUBSCRIPTS 

2 

w 

CD 

i.  j 


Conditions  immediately  downstream  of  normal  shock  wave 

Sphere  surface 

Free -stream  conditions 

Shadowgraph  station  interval,  i  =  1  to  6,  j  =  2  to  6 


SECTION  i 
INTRODUCTION 
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Falling,  spherically  shaped,  pressurized  balloons  are  extensively 
used  for  atmospheric  probing  (Ref,  1).  These  probes  can  be  used  to 
determine  several  properties  of  the  upper  atmosphere  (horizontal  winds, 
density,  temperature,  and  pressure),  the  most  important  property  being 
atmospheric  density.  The  velocity  of  descent  of  a  sphere  of  specified 
size  and  weight  at  a  particular  altitude  is  a  function  of  the  density  of  the 
atmosphere  at  that  altitude  and  the  drag  coefficient  of  the  sphere  at 
those  conditions.  Density  is  inversely  >roportional  to  the  balloon  drag 
coefficient  and  can  be  inferred  if  that  drag  coefficient  is  known. 

For  several  years  the  Air  Force  Cambridge  Research  Laboratory 
has  conducted  falling  sphere  atmospheric  soundings  between  the  altitudes 
of  140  and  30  km.  Two  different  kinds  of  balloons  have  been  used,  and 
the  motion  of  each  is  measured  by  a  different  method.  One  probe  is  a 
1-m-diam  Mylar® passive  sphere,  ROBIN,  whose  space-time  positions 
during  its  fall  are  measured  with  precision  ground-based  radar.  The 
second  balloon  is  an  "Instrumented  Sphere"  whose  falling  motion  is 
measured  by  on-board  accelerometers.  These  two  balloons  have  dif¬ 
ferent  fall  rates  and  trajectories.  Figure  1  (Appendix  I)  is  a  Mach 
number-Reynolds  number  envelope  of  typical  trajectories  of  the  two 
types  of  balloons. 

The  need  for  further  sphere  drag  measurements  arose  from  two 
sources,  namely: 

1.  The  results  of  previous  sphere  drag  tests,  by  many 
experimenters  and  using  a  variety  of  test  facilities,  are 
not  always  in  agreement. 

2.  Falling  spheres  experience  some  M^-Re,,,  values  for 
which  drag  coefficients  are  not  available. 

It  is  fo»~  the  reasons  stated  above  that  the  present  sphere  drag  study  was 
undertaken.  The  investigation  has  been  conducted  in  the  M^-Re^  range 
defined  in  Fig.  1  according  to  the  AFCRL  specified  priorities. 

The  present  sphere  drag  investigation  was  undertaken  in  the 
von  Karman  Gas  Dynamics  Facility  (VKF)  of  AEDC.  In  an  aeroballistic 
range,  sphere  drag  coefficient  was  measured  for  a  broad  range  of  Mach 
number  (0.  12  <  M,,,  <  6.  39)  and  Reynolds  number  (15  <  Re^  <  50,  300), 
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SECTION  II 

AEROBALLISTIC  RANGE 


The  aeroballistic  range  testing  was  conducted  in  the  Hyperballistic 
Range  (K).  This  range  is  a  variable  density,  free -flight  test  unit  that 
is  used  for  aerophysical  testing  and  for  classical  aerodynamic  tests, 
notably  the  measurement  of  drag.  The  aeroballistic  range  consists 
basically  of  two  tanks,  a  pumping  system,  a  launcher,  and  several  in¬ 
strumentation  systems. 


2.1  RANGE  K 

The  range  tank  and  blast  tank  are  6-ft-diam.  steel  cylinders  con¬ 
nected  by  a  short  spool  piece  containing  a  high  vacuum  valve  which  per¬ 
mits  isolation  of  the  two  tanks  (see  Fig.  2).  The  14-ft-long  blast  tank 
is  situated  between  the  launcher  and  he  range  tank.  The  blast  tank  is 
used  to  absorb  the  expanding  muzzle  gases,  and  it  is  in  this  chamber 
that  the  model  is  separated  from  the  sabot  which  adapts  it  to  the  bore 
of  the  launcher.  Once  separated  from  the  sabot,  which  is  stopped  in 
the  blast  tank,  the  model  flies  through  a  3-in.  -diam  hole  in  the  spool 
piece  and  into  the  range  tank. 

The  range  tank  is  100  ft  long  and  is  equipped  with  the  instrumenta¬ 
tion  necessary  to  make  the  desired  test  measurements.  Shadowgraph, 
schlieren,  temperature,  and  pressure  instrumentation  systems  are 
permanently  installed.  Flight  is  terminated  by  model  impact  on  a  thick 
plate  at  the  end  of  the  range  tank. 

A  three-stage  vacuum  pumping  system  provides  the  desired  range 
pressure.  The  blast  and  range  tanks  have  independent  pumping  sys¬ 
tems  which  facilitate  testing  at  low  pressures  because  the  range  tank 
can  be  kept  isolated  at  low  pressure  while  the  launcher  is  being  pre  ¬ 
pared  for  the  next  shot.  For  this  test,  all  shots  were  fired  with  the 
blast  and  range  tank  pressures  equal. 

Three  different  launchers  were  employed  in  achieving  the  desired 
range  of  velocities,  0.  2  <  M,,,  <6.0: 

1.  Cold-Gas  Gun  -  Most  of  the  spheres  for  which  <1.7 
were  launched  with  a  single-stage,  cold-gas  pneumatic 
launcher  using  helium  as  the  driving  gas.  This  launcher 
is  essentially  the  same  as  that  used  by  Lawrence  (Ref.  7) 
in  an  earlier  VKF  sphere  drag  test. 
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2.  2o-mm  Cannon  -  The  majority  of  the  spheres  for  which 

>1.7  were  launched  with  a  20-mm-diam,  smooth¬ 
bore,  single-stage  launcher.  Burning  gunpowder  sup¬ 
plied  the  driving  pressure. 

3.  Rule  -  A  few  shots  were  made  with  a  30-cal,  rifled- 
bore,  single-stage,  powder-driven  launcher. 

Launchers  1  and  2  were  used  for  92  percent  of  the  shots. 


2.2  RANGE  INSTRUMENTATION 

Six  dual-axis  spark  shadowgraph  stations  record  the  model  position¬ 
time  history  for  each  flight.  These  stations  are  located  at  nominally 
1 5  -  ft  intervals  along  the  range  tank.  Optical  axes  at  each  station  are 
mutually  orthogonal  with  the  range  centerline.  The  shadowgraph  con¬ 
sists  of  a  pinhole  spark  source  on  one  side  of  the  range  and  a  Fresnel 
lens  and  camera  on  the  other  side.  Except  for  the  Fresnel  lens,  the 
shadowgraph  system  components  are  external  to  the  range  tank.  The 
spark  has  a  duration  of  0.  1  /usee;  therefore,  this  is  the  effective  ex- 
posuie  time  of  the  shadowgram.  A  shadow  detector  is  used  to  trigger 
the  spark  when  the  model  is  near  the  center  of  the  shadowgraph  field  of 
view.  When  the  spark  is  initiated,  a  signal  is  sent  to  a  multichannel 
digital  event  chronograph  which  iccords  the  times  of  successive  shadow- 
gram  exposure.  The  Fresnel  lenses  bear  scribed  fiducial  grids,  and 
the  stations  have  beer,  surveyed  accurately  as  to  orthogonality  and  spacing 
so  that  model  position  can  be  determined  in  real  space  with  respect  to 
master  range  axes.  This  survey  has  been  further  refined  by  calibration 
shots  with  spheres  of  high  density  material  fired  into  a  hard  vacuum 
(AV  =  0)  at  velocities  which  minimize  time -position  uncertainties. 

Gas  temperatures  in  the  range  tank  are  measured  by  means  of 
mercury-in-glass  thermometers  and  by  copper-constantan  thermo¬ 
couples  read  out  on  a  strip-chart  recorder,  A  pair  of  these  probes  is 
located  near  each  of  the  six  shadowgraph  stations. 

In  order  to  ensure  that  the  air  in  the  range  is  sufficiently  dry,  dew¬ 
point  temperature  is  measured  with  a  hygrometer  located  near  the  center 
of  the  range  tank.  Dry  air  can  be  supplied  by  the  VKF  high  pressure 
bottle  (H2O  <  100  ppm)  or  from  a  source  of  pure  breathing  air 
(H2O  <  30  ppm). 

A  high  sensitivity,  single-pass  schlieren  photographic  system  is 
available  for  visualizing  the  model  flow  field.  The  viewfield  diameter 
is  12  in.,  and  the  system  is  operable  in  either  single-frame  or  multiple- 
frame  modes. 
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All  gages  used  to  measure  pressure  in  the  range  are  connected  to  a 
large  stainless  steel  manifold  which  can  be  isolated  from  the  range  tank 
and  kept  under  hard  vacuum  when  not  in  use.  As  many  as  six  gages  can 
be  connected  to  the  manifold.  For  this  test  three  types  of  gages  were 
used: 

1.  Hass  Mercury  Manometer  -  This  instrument  is  used 
over  the  pressure  range  from  15  to  760  torr.  It  has 
a  resolution  of  0.  05  torr. 

2.  Micromanometer  -  This  is  an  oil-filled,  U-tube 
manometer  with  which  VKF  has  considerable  experi¬ 
ence.  It  is  used  for  measuring  pressures  between  0.  1 
and  15  torr.  Resolution  of  the  micromanometer  is 

0. 00075  torr. 

3.  Baratron®  -  This  is  a  variable-capacitance  pressure 
transducer  which  has  a  small  internal  volume.  The 
volume  of  the  other  two  instruments  is  large  and  they, 
therefore,  do  not  respond  as  rapidly  as  the  Baratron  to 

a  changing  pressure.  The  Baratron  has  two  sensing  heads: 

1  torr  and  30  torr.  This  permits  the  measurement  of 
pressures  from  30  torr  down  to  approximately  10-^  torr. 
Baratron  resolution  is  1  x  10"^  torr  for  pressures  above 
10  torr  and  is  1  x  10”^  torr  for  pressures  below  10  torr. 

Insofar  as  possible,  these  pressure  gages  are  calibrated  with  standards 
traceable  to  the  National  Bureau  of  Standards  (NBS).  In  pressure  regions 
where  no  NBS  standard  exists  (p  <  1  torr),  calibration  is  accomplished 
by  reasonable  extension  of  existing  standards  and  by  means  of  a  VKF- 
built  calibrator  described  in  Ref.  3.  Since  Ref.  3  was  written,  a  new 
calibrator  has  been  built  to  closer  tolerances,  and  greater  precision  is 
now  possible.  The  calibr  ator  itself  has  been  thoroughly  checked  against 
the  NBS  standards  where  possible.  Calibration  of  all  the  gages  was  per¬ 
formed  before,  during,  and  after  this  test  in  order  to  ensure  the  greatest 
possible  accuracy  in  the  pressure  measurement. 

2.3  MODELS  AND  SABOTS 

The  variety  of  sphere  mater  ials  and  sizes  used  in  the  range  test  are 
listed  in  Table  I  (Appendix  II)  and  illustrated  in  Fig.  3.  Whenever  pos¬ 
sible,  spheres  of  a  conventional  material  were  used.  Commercially 
available,  precision  grade,  standard  spheres  were  chosen  in  order  to 
ensure  high  quality  surface  finish  and  to  minimize  nonsphericity.  The 
spheres  were  then  individually  inspected  and  measured  at  AEDC.  The 
variety  of  spheres  used  was  necessary  in  order  to  produce  the  t  -  aired 
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Hcvnohls  number  aiul  Mach  number  at  a  velocity  drop,  AV,  of  between  1 
and  (i  percent  of  the  average  velcu-ity,  V. 

I'lie  minimum  model  size  of  3/32-in,  diameter  was  the  smallest 
that  could  he  i  eliably  detected  by  the  present  shadowgraph  triggers. 

This  size  limitation,  along  with  the  desire  to  maintain  AV  between  1 
and  ti  percent,  establishes  the  minimum  Reynolds  number  at  which  tests 
can  lie  conducted  with  conventional  model  materials.  Models  of  molded 
Dylite®  were  used  to  overcome  this  limitation.  Dylite  is  a  foamed 
plastic  that  was  used  by  VKF  in  the  manufacturing  of  spheres  whose 
density  ranged  from  1  to  (i  lb/ft**.  References  4  and  5  discuss  the 
original  use  of  Dylite  in  the  manufacturing  of  lightweight  models  and  the 
measurement  of  C[)  with  these  models.  The  use  of  these  lightweight 
foam  models  allowed  testing  at  lower  Re,,,  while  maintaining  the  desired 
model  deceleration.  The  lower  Re^  limitation  for  Dylite  models  was  not 
reached  in  this  test. 

Several  sabot  designs  were  used  during  the  course  of  this  test. 
Typical  examples  are  shown  in  Fig.  3.  Except  for  a  few  Dylite  and 
partially  Dylite  sabots,  all  were  made  of  Lexan®.  A  different  method 
of  separating  the  sphere  from  the  sabot  was  employed  with  each  of  the 
three  launchers.  Figure  4  illustrates  these  stripping  methods. 


SECTION  III 

RANGE  SPHERE  DRAG  MEASUREMENT 


3.1  THEORY 


The  equation  of  motion  of  a  sphere  in  free  flight  is  simply  Newton's 
Law  written  as  follows: 


D  =  ma  =  1/2  p  V2CDA 

Using  the  relationships  p  =  pRT  and  A  =  7r/4  d2  and  substituting  in 
Eq.  (1),  it  can  now  be  written 

dV  *  P  V  CD  d 
m  dt  =  8  RT 


(1) 


(2) 


or,  solving  for  Cjy 


'D 


-8™  3  (dV/dt) 

7T  p  d2  V2 


(3) 


Making  the  substitution 

1/V  (dV/dt)  =  dt/dx  •  dV/dt  =  dV/dx 
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Eq.  (3)  takes  the  form 

CD  =  (dV/dx)  (4) 

it  p  V 

It  is  shown  that  drag  coefficient  can  be  obtained  if  the  sphere  mass, 
diameter,  velocity,  and  deceleration  and  the  range  temperature  and 
pressure  are  known.  These  quantities  can  all  be  measured  in  a  range 
experiment. 


3.2  METHOD 

Before  launching,  the  sphere  is  cleaned,  inspected,  weighed,  and 
measured  to  determine  its  surface  finish  quality,  m?”  and  diameter. 
The  sphere  is  then  launched  into  dry  air  at  a  known  p.  assure  and  tem¬ 
perature.  Position-time  measurements  are  obtained  during  the  flight 
of  the  sphere,  and  then  velocity  and  deceleration  are  derived  from  these 
data.  The  sphere  mass,  diameter,  and  material  and  the  range  pressure 
are  varied  in  order  to  produce  the  desired  Reynolds  number  at  a  given 
velocity  and  to  keep  the  velocity  drop  between  1  and  6  percent. 


3.3  DATA  REDUCTION 

The  sphere  image  on  the  shadowgrams  is  read  with  a  digitized  film 
reader,  and  the  resultant  numerical  values  are  converted.,  by  means  of 
a  computer  program  containing  the  range  shadowgraph  system  survey, 
into  real  space  positions.  These  model  positions  and  the  associated 
timing  values  are  used  to  calculate  the  average  velocity  between  all  pairs 
of  adjacent  shadowgraphs  (Vj^i+i  for  i  =  1  to  5).  In  addition,  average 

velocities  over  a  two-station  interval  (Vj  i+2  f°r  i  =  1  to  4)  are  calcu¬ 
lated.  This  second  sot  of  velocities  has  only  half  of  the  uncertainty  in  x 
and  t,  compared  to  the  first  set,  since  the  total  distance  and  time  are 
doubled. 

A  least-squares  straight  line  is  computer  fitted  to  the  velocity- 
position  data.  Figure  5  is  a  typical  example  of  this  fitting.  By  fitting  a 
straight  line  to  these  data,  the  assumption  is  made  that  Cp>  is  a  constant 
during  the  flight.  This  assumption  is  quite  valid  for  this  test  since  the 
model  velocity  drop,  AV,  was  less  than  6  percent  of  the  midrange 
velocity.  The  value  of  V  used  in  the  computation  of  Re^  and  Cp>  and  the 
value  of  dV/dx  used  in  the  computation  of  Cp  are  taken  from  the  fitted 
line.  In  this  data  reduction,  Eq.  (4)  is  used  to  calculate  Cp). 
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Another  data  reduction  method  is  also  used  on  the  position-time 
data.  The  time-position  relationship  of  a  model  in  free  flight  can  be 
represented  by  a  cubic  equation  in  distance,  viz 

t  =  aQ  -raj  x  +  a;>  +  a3  x^ 

This  equation  is  computer  fitted  to  the  time -position  data  by  the  method 
of  least  squares.  In  this  case  the  deceleration  is  not  assumed  constant. 

If  the  model  undergoes  a  large  percentage  velocity  drop,  this  treatment 
of  data  is  required.  Reynolds  number  and  Cq  are  computed  from  values 
of  V  and  dV/dt  taken  from  the  fitted  cubic.  In  this  treatment,  Cg>  is 
calculated  from  Eq.  (3). 

Although  the  cubic  fit  was  not  necessary  for  the  treatment  of 
position-time  data  on  this  test  since  0.  01  <  AV/V  <  0.  06,  both  linear  fit 
derived  and  cubic  fit  derived  values  of  M^,  Re^,  and  Cg)  were  computed 
for  all  shots.  To  verify  the  contention  that  the  linear  fit  is  adequate  if 
AV/V  <  0.  06,  differences  in  linear  and  cubic  fit  derived  values  of 
Mm,  Re^,  and  Cg)  were  obtained  for  seven  shots  where  0.058  <  AV/V  < 

0.  073.  These  shots  were  picked  because  these  differences  will  increase 
with  increasing  AV/V.  For  all  seven  shots,  the  differences  in  and  Rero 
were  less  than  0.  01  percent  and  for  Cd  were  less  than  0.  1  percent. 

3.4  MEASUREMENTS  AND  UNCERTAINTY 

3.4.1  Distance 

As  explained  in  Section  4,  3,  downrange  position  is  determined  from 
the  shadowgrams.  The  results  of  past  Range  K  shadowgraph  station  sur¬ 
veys  and  calibration  shots  in  conjunction  with  check  calibration  shots 
fired  during  the  course  of  this  test  show  that  the  maximum  error  in  model 
position  (flight  distance)  is  less  than  ±0.  0025  in  The  flight  distance  be¬ 
tween  adjacent  shadowgraph  stations  is  nominally  15  ft. 

3.4.2  Time 

The  digital  event  chronograph  associated  with  the  shadowgraph  sys¬ 
tem  provides  timing  values  with  a  resolution  of  1  x  10”^  sec.  The  maxi¬ 
mum  error  in  the  time  of  flight  between  adjacent  shadowgraphs  is 
±1  x  10'^  sec  (±0.  1  m sec).  For  this  test,  interstation  flight  time  varied 
between  1.  1  x  10^  Msec  (M^  =  0.  12)  and  2.  1  x  10^  m sec  (MB  =  6.  39).  The 
timing  error  is  not  cumulative.  The  maximum  error  in  the  total  time  of 
flight  is  also  ±0.  1  Msec. 
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3.4.3  Velocity 

The  velocity  used  in  the  computation  of  Cp>  and  Re^  is  the  midrange 
velocity  as  derived  from  the  linear  fit  to  the  Vi#  j  and  x  data.  The  maxi¬ 
mum  error  in  is  0.03  percent.  On  this  test,  shots  were  rejected 

if  the  maximum  computed  from  x,  t  data)  -  from  curve 

fit)]  was  more  than  0.  025  percent  of  the  midrange  velocity.  This  maxi¬ 
mum  velocity  error  was  less  than  <0.  02  percent  on  95  percent  of  the 
shots.  The  midrange  velocity,  used  in  calculating  Cjy  and  Pe^,  is  known 
to  an  uncertainty  of  less  than  0.  02  percent,  since  it  is  obtained  from  a 
curve  fit  to  all  the  velocity-position  data.  The  velocity  range  on  this 
test  was  from  135  to  7250  ft/sec. 

3.4.4  Deceleration 

Deceleration  is  one  of  the  least  accurately  known  quantities.  Where¬ 
as  small  perturbations  around  the  true  values  of  x  and  t  produce  very 
small  changes  in  velocity,  such  perturbations  can  produce  a  relatively 
large  change  in  dV/dx.  An  analysis  of  x,  t,  and  dV/dx  errors  indicates 
that  for  the  v/ide  range  of  shots  made  during  this  test,  the  maximum 
possible  uncertainty  in  deceleration  can  be  greater  than  1  percent. 
However,  the  x,  t,  and  V  accuracy  on  the  majority  of  the  shots  is  such 
that  the  uncertainty  in  dV/dx  should  not  exceed  1  percent. 

3.4.5  Mass 

Sphere  mass  was  determined  on  either  of  two  analytical  balances  or 
on  a  torsion  balance.  The  choice  was  dependent  upon  sphere  mass. 

Great  care  was  taken  in  the  handling  of  the  spheres,  and  they  were 
weighed  after  they  had  been  cleaned  and  their  diameters  measured. 
Balance  zero  was  checked  before  each  weighing,  and  each  sphere  was 
weighed  several  times.  For  spheres  of  m  <  0.  021  gm  (those  weighed  on 
the  torsion  balance),  a  class  M  standard  weight  of  approximately  equal 
mass  was  weighed  before  and  after  each  sphere  weighing.  Uncertainties 
in  the  measurement  of  sphere  mass  are  as  follows: 

0.  00023  gm  for  m  >  0.  5  gm 
0.  00016  gm  for  0.  14  <  m  <  0.  5  gm 
0.  00005  gm  for  0.  05  <  m  <  0.  14  gm 
0.  00002  gm  for  0.  015  <  m  <  0.  05  gm 
0.  00001  gm  for  m  <  0.  015  gm 


3.4.6  Diameter 

Two  instruments  were  employed  in  the  measurement  of  sphere 
diameter.  Spheres  of  conventional  materials  were  measured  with  a 
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light  wave  micrometer  which  has  an  accuracy  of  ±0.  00001  in.  Since 
this  instrument  applies  mechanical  pressure  to  the  sphere  during 
measurement,  this  technique  is  not  suitable  for  the  foam  spheres. 

These  were  measured  with  a  microscope  which  has  an  accuracy  of 
±0.  0001  in.  Several  measurements  were  made  on  each  sphere  in  order 
to  determine  the  maximum  and  the  minimum  diameters.  The  diameter 
used  in  the  calculation  of  Cd  and  Re^,  was  the  average  of  these  two 
diameters.  Spheres  of  conventional  materials  were  selected  such  that 
the  uncertainty  in  diameter  attributable  to  nonsphericity  was  always 
less  than  0.  1  percent  of  the  diameter.  For  the  foam  spheres,  the  un¬ 
certainty  in  diameter  attributable  to  nonsphericity  was  0.  5  percent  for 
d  =  0.  25  in.  and  1.  0  percent  for  d  =  0.  125  in. 

3.4.7  Temperature 

Range  temperature,  at  any  probe  location,  is  most  accurately 
determined  by  the  mercury-in-glass  thermometer  which  has  a  maxi¬ 
mum  error  of  ±0.  1°F.  The  thermocouple  temperatures  have  an 
accuracy  of  ±0.  5°F  and  were  used  as  a  check  on  the  thermometer  tem¬ 
perature  readings  in  order  to  determine  that  there  were  no  transients 
at  the  time  of  the  shot.  Even  without  noticeable  transients  at  any  pr  _>be 
location,  there  usually  exists  a  slight  gradient  in  temperature  along  the 
length  of  the  range.  The  temperature  at  any  probe  location  never  dif¬ 
fered  from  the  average  of  the  six  temperature  measurements  by  more 
than  ±2°F.  This  maximum  variation  yields  an  uncertainty  of  no  greater 
than  0.  38  percent  in  absolute  temperature.  Average  range  temperature 
on  this  test  varied  between  extremes  of  6  3.  2  and  83.  8°  F,  but  was 
normally  between  73  and  79°  F. 

3.4.8  Pressure 

The  pressure  variation  of  this  test  was  from  5  38  to  0.  034  torr. 

This  required  the  use  of  all  three  types  of  gages  listed  in  Section  2.  2. 
The  use  of  each  gage  was  limited  to  that  pressure  range  over  which  it 
had  the  greatest  accuracy.  Where  pressure  ranges  for  more  than  one 
gage  overlapped,  each  was  read.  The  difference  in  pressure  between 
gages  was  no  greater  than  1. 0  percent  even  at  the  end  of  a  particular 
gage's  range.  This  does  not  necessarily  imply  a  pressure  uncertainty 
of  1.0  percent  since  the  pressure  reading  of  the  gage  nol  at  the  extreme 
of  its  pressure  range  would  be  the  more  accurate.  The  maximum  un¬ 
certainty  in  pressure  ir  no  greater  than  1  percent  for  p  >  1  torr.  For 
pressures  below  1  torr,  the  uncertainty  increases  with  decr  easing 
pressure  and  reaches  a  maximum  no  greater  than  2  percent  at 
p  -  0.  034  torr. 
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3.4.9  Drag  Coefficient 

The  uncertainties  specified  above  are  maximum  possible  uncer¬ 
tainties.  Maximum  probable  uncertainties  would  be  less  than  or  equal 
to  these  values.  The  maximum  probable  uncertainty  in  drag  coefficient 
attributable  to  uncertainties  in  the  measurements  are  as  follows: 

V  t  02  percent 

T  0. 28  percent 

dV  /dx  1.0  percent 

m  0.  1  to  0.  3  percent  (conventional  materials) 

0.5  percent  ( 1 /4-in. -diam  Dylite) 

1.  0  percent  ( 1  /  8-in.  -diam  Dylite) 

d  0.  2  to  0.  3  percent  (conventional  materials) 

1.  0  percent  ( 1/4-in.  -diam  Dylite) 

2.  0  percent  ( 1/ 8-in.  -diam  Dylite) 

p  0.  5  percent  (p  >1  torr) 

1.  0  percent  (0.  1  <  p  <  1  torr) 

2.  0  percent,  (p  <  0.  1  torr) 

As  can  be  seen  from  the  data  in  Figs.  11a  through  h  and  111  through  v, 
the  scatter  around  the  fitted  curve  is  such  that  Cp)  is  certainly  deter¬ 
mined  to  within  ±2  percent  for  Re,,,  >  200.  Below  Re,,,  =  200,  the  scatter 
increases  because  of  increasing  pressure  uncertainty  for  all  conditions, 
but  the  scatter  becomes  noticeably  worse  only  when  Dylite  models  are 
used.  In  spite  of  the  increased  data  scatter,  the  greater  number  of 
shots  fired  at  these  conditions  and  the  method  of  data  treatment 
(Section  5.  1)  results  in  values  of  Cp)  that  are  probably  correct  within 
±2  percent  even  down  to  Re*  -  20. 

3.4.10  Reynolds  Number 

The  accuracy  of  Re*  is  dependent  upon  the  uncertainties  in  p  ,  T*,  V, 

QD 

and  d.  The  maximum  probable  uncertainties  in  Re*  are  as  follows: 

0.  8  percent  for  Re*  >  200 
1.  2  percent  for  Re*  <  200 


3.4.11  Mach  Number 

As  M*  depends  only  upon  V  and  T,  the  maximum  probable  uncer¬ 
tainty  is  0.  2  percent. 
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4.1  AEROBALLISTIC  RANGE  DATA 

The  results  of  the  present  test  program  for  0.  1  <  Mm  <6.0  and 
2  x  10*  <  Re*  <  105  are  listed  in  Table  II  and  are  presented  in  Figs.  6a 
through  v  as  the  variation  of  drag  coefficient  with  free-stream  Reynolds 
numbers  for  discrete  Mach  number  intervals.  Also  shown  in  these 
figures  are  the  results  of  sphere  drag  measurements  presented  in 
Refs.  2,  5,  and  6  through  12.  No  attempt  has  been  made  to  ensure  that 
this  listing  of  references  is  complete.  These  data  were  selected  be¬ 
cause  the  ratio  of  the  model  wall  temperature  to  free-stream  tempera¬ 
ture  was  approximately  unity,  as  is  the  case  for  the  ballistic  range  data 
contained  herein. 

At  subsonic  speeds  (Figs.  6a  through  j),  the  present  data  are  in 
good  agreement  with  those  obtained  by  Lawrence  (Ref.  2)  in  the  same 
facility.  The  data  contained  in  Refs.  6  and  7  for  subsonic  velocities  are 
consistent  with  the  data  obtained  in  the  present  investigation  (Figs.  6a 
through  j).  For  supersonic  velocities,  the  data  available  from  a  variety 
of  sources  (Refs.  5  and  6  through  11)  are  also  consistent  with  those  ob¬ 
tained  in  the  present  investigation  (Figs.  6k  through  v). 

Curves  have  been  faired  through  all  of  the  available  data  at  the 
discrete  Mach  number  intervals  listed  in  Figs.  6a  through  v.  An  aver¬ 
age  Mach  number  has  been  assigned  to  each  of  these  Mach  number 

intervals  (e.  g.  ,  for  interval  0.  19  <  Mm  <  0, 27,  Mm  =  0.  23)  and  cross- 

®avg 

plots  of  Cp)  versus  Mach  number  have  been  derived  at  fixed  Reynolds 
numbers  (Figs.  7a  through  d).  The  symbols  shown  in  Fig.  7  are  not 
experimental  points  but  are  values  derived  from  faired  curves  in  Figs.  6a 
through  v.  It  is  from  these  crossplots  that  drag  coefficients  should  oe 
derived.  An  estimation  of  Cp>  at  -  0  is  shown  in  Appendix  IV. 

It  is  evident  from  a  consideration  of  the  data  contained  in  Figs.  6a 
through  v  that  for  Re,,,,  <  100  the  spread  in  the  experimental  values  of  Cp> 
increases.  In  general,  measurements  in  this  regime  were  obtained  with 
ultralight  spheres  and  at  ambient  pressures  less  than  100  pHg.  Several 
factors  contribute  to  this  spread  --  larger  uncertainties  in  model  weight 
and  diameter,  model  distortion,  and  pressure. 

From  an  examination  of  the  measurement  errors  discussed  in 
Section  4.4  and  the  consistency  of  the  experimental  data  (Figs.  6a 
through  v),  indications  are  that  the  total  errors  in  sphere  drag  coeffi¬ 
cient  derived  from  the  faired  curves  are  no  greater  than  ±2  percent. 
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4.2  SPHERE  DRAG  COEFFICIENTS  AT  SUBSONIC  SPEEDS 

An  accurate  knowledge  of  the  drag  coefficient  of  a  sphere  at  low 
subsonic  speeds  is  of  interest  in  the  ROBIN  Falling  Sphere  program. 
From  a  consideration  of  the  possible  falling  sphere  trajectories  (Fig.  1), 
a  knowledge  of  the  sphere  drag  coefficients  for  0.  1  <  M,,,  <  0.  5  and 
103  <  Re*  <  5  x  104  is  required. 

Before  1930,  many  measurements  of  the  drag  of  a  sphere  falling 
through  various  fluids  were  made,  and  a  body  of  information  was 
generated  for  10“ *  <  Re*,  <  103.  It  was  considered  that  these  data  gave 
the  level  of  the  incompressible  drag  coefficient  of  a  sphere  in  steady 
nonturbulent  flow.  Usually  these  data,  which  have  a  significant  degree 
of  scatter,  are  represented  by  a  single  line  which  is  called  the  "standard 
drag  curve.  "  This  curve  has  appeared  in  many  textbooks  and  reports  at 
least  as  far  back  as  1931  (Ref.  13).  (No  attempt  is  made  here  to  deter¬ 
mine  its  exact  origin. ) 

The  "standard  drag  curve"  for  the  range  of  interest  of  the  ROBIN 
Spheres  is  shown  in  Fig.  8.  Also  shown  in  this  figure  are  two  examples 
of  the  early  data  upon  which  this  curve  has  been  based  (Refs.  14  and  15). 

Heinrich,  Niccum,  and  Haak  (Ref.  16)  made  some  sphere  drag 
measurements  in  a  wind  tunnel  for  0.  078  <  M,,,  <  0.  39  and  2  x  103  < 

Re*  <  2  x  104.  These  data  are  compared  with  the  standard  drag  curve 
in  Fig.  8  and  are  shown  to  be  significantly  higher  than  the  standard 
values.  Sivier  (Ref.  17)  has  measured  the  drag  of  magnetically  sup¬ 
ported  spheres  in  a  wind  tunnel  with  a  free -stream  turbulence  intensity 
up  to  b  percent.  These  sphere  drag  values  (cf.  Fig.  8)  are  also  signifi¬ 
cantly  greater  than  the  standard  drag  values.  Zarin  (Ref.  18)  refined 
the  magnetic  balance  system  used  by  Sivier  (Ref.  17)  and  varied  the 
free-stream  turbulent  intensity  level.  He  obtained  some  sphere  drag 
measurements  at  free-stream  turbulent  intensity  levels  of  less  than 
1  percent.  These  values  are  shown  in  Fig.  8.  For  Re^  >  103  these 
values  are  still  significantly  greater  than  the  standard  values.  For 
Re®  <  103  these  values  are  in  good  agreement  with  the  standard  values. 
From  his  study,  Zarin  concluded  that  small  degrees  of  free-stream 
turbulence  were  the  cause  of  his  higher -than-standari  drag  values. 

Also  shown  in  Fig.  8  are  the  present  drag  values  from  Fig.  fib.  It 
can  be  seen  that  these  values  are  in  reasonable  agreement,  with  the 
standard  values  for  5  x  10^  <  Re„,  <  104.  It  is  reasonable  to  assume 
that  the  ballistic  range  data  are  representative  of  a  free-stream  turbu¬ 
lent  intensity  level  approaching  zero.  If  the  effects  of  turbulence  shown 
by  Zarin  Ref.  18)  are  correct,  then  it  might  be  inferred  that  some 
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consideration  of  the  scaling  of  free -stream  turbulence  intensity  level 
may  be  necessary  before  sphere  drag  coefficients  from  turbulent  test 
facilities  can  be  applied  to  a  particular  falling  sphere  in  the  atmosphere. 
It  is  an  experimentally  demonstratable  fact  that  turbulence  can  affect 
the  critical  Reynolds  number  where  boundary-layer  transition  occurs. 
The  critical  Reynolds  number  is  that  Reynolds  number  for  which  the 
drag  coefficient  decreases  sharply  and  passes  through  Cp>  =  0.  3  as 
Re*  increases.  Increasing  degrees  of  turbulence  have  been  shown  to 
reduce  the  critical  Reynolds  number  for  Re,.  >10^  (Rc-fs.  19  and  20) 
and  to  increase  Cp)  in  the  area  100  <  Re*  <  1000  (Ref.  18).  Millikan 
and  Klein  (Ref.  19)  found  that  the  critical  Reynolds  number  for  a  sphere 
in  free  flight  in  the  atmosphere  was  not  affected  by  varying  atmospheric 
turbulence  structure.  This  result  was  explained  by  the  fact  that  the 
scale  of  the  turbulence  in  the  atmosphere  is  large  compared  with  the 
boundary  layer  on  the  sphere.  It  would  seem  reasonable  to  assume  that 
if  the  critical  Reynolds  number  of  a  free -flight  sphere  is  unaffected  by 
the  full-scale  turbulence,  then  the  drag  coefficient  at  subcritical 
Reynolds  numbers  would  also  be  unaffected  by  turbulence. 

As  a  matter  of  interest,  there  appear  to  be  inconsistencies  in  the 
literature  as  to  the  effect  of  free-stream  turbulence  upon  drag  coeffi¬ 
cient.  Zarin  (Ref.  18)  indicates  that  his  and  Sivier's  (Ref.  17)  greater 
than  "standard  values"  for  sphere  drag  can  be  explained  in  terms  of 
turbulence,  intensity,  and  scale.  However,  Probstein  and  Fassio 
(Ref.  21)  explain  Ingebo's  (Ref.  22)  lower  than  standard  values  of  drag 
in  terms  of  turbulence  in  the  flow.  Actually,  depending  on  the  particular 
Reynolds  number  regime,  it  would  seem  that  turbulence  could  either 
raise  or  lower  Cp). 

4.3  EFFECTS  OF  COMPRESSIBILITY  ON  SPHERE  DRAG  COEFFICIENT 

As  in  the  case  for  the  subsonic  velocities,  there  is  a  generally 
accepted  curve  (Fig.  9)  which  has  been  used  to  indicate  the  effects  of 
compressibility  upon  sphere  drag  coefficient.  This  curve  appears  to 
have  originated  with  Hoerner  (Ref,  23)  and  has  been  repeated  in  several 
references  since  (e.  g. ,  Refs.  1,  17,  and  18).  The  form  of  the  varia¬ 
tion  of  sphere  drag  coefficient  with  Mach  number  (Ref.  23)  was  based 
upon  several  sets  of  data  for  0  <  M^  <1.0  obtained  before  1946  (in¬ 
cluding,  for  example,  the  ballistic  range  data  of  Ref.  6).  The  curve  is 
characterized  by  a  pronounced  dip  in  the  drag  coefficient  at  -  0.  85. 
Unfortunately,  the  ballistic  range  test  (Ref,  6)  did  not  obtain  data  for 
0.  65  <  M*  <  0.  85.  An  inspection  of  the  remaining  data  (Ref.  23)  upon 
which  this  curve  (Fig.  9)  was  based  indicates  that  the  curve  fitted  to  the 
data  may  not,  in  fact,  be  the  best  fit  to  the  data.  This  is  confirmed  to 
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some  extent  when  Naumann's  data  (Ref.  7)  is  considered  in  Fig.  9.  It 
can  be  seen  (Fig.  9)  that  these  data  agree  with  Charters'  and  Thomas' 
data  (Ref.  6)  and  do  not  show  a  dip  in  drag  coefficient  at  MB  -  0.  85. 

The  reason  for  presenting  the  above  discussion  is  to  examine  the 
validity  of  the  dip  in  light  of  the  present  results.  If  the  form  of 
Hoerner's  curve  is  accepted,  then  to  some  extent  his  explanation  for 
the  reason  for  the  dip  is  accepted  (cf.  Ref.  18).  Hoerner's  explanation 
for  the  dip  is  that  it  is  attributable  to  "a  favorable  interaction  between 
the  local  supersonic  field  of  flow  existing  at  and  behind  the  location  of 
the  cylinder's  (sphere's)  maximum  thickness  and  the  flow  pattern  within 
its  wake.  "  If  this  reasoning  is  accepted,  then  for  lower  subcritical 
Reynolds  numbers  a  curve  with  similar  characteristics  would  be  ex¬ 
pected.  None  of  the  summary  curves  shown  in  Fig.  7  indicates  that 
this  occurs.  To  indicate  that  the  effect  has  not  been  obscured  by  the 
smoothing  procedure  used  in  deriving  the  summary  curves,  the  data 
obtained  for  1930  <  Re^  <  2080  are  shown  in  Fig.  10.  These  and 
Naurnann's  (Ref.  7)  data  are  sufficiently  well  defined  to  state  with  some 
certainty  that  there  is  no  dip  in  the  Cq  versus  M,,,  curve  for  Mffl  «  0.  85. 

4.4  TRANSONIC  SPHERE  DRAG  COEFFICIENTS 

There  appear  to  be  very  few  experimental  measurements  of 
sphere  drag  in  the  transonic  speed  regime.  The  data  contained  in 
Figs.  6a  through  1  and  repeated  in  the  summary  curves  in  Figs.  7a 
through  d  represent  the  most  comprehensive  set  of  results  in  this 
speed  regime.  From  a  consideration  of  the  data  presented  in  Fig.  7 
for  2  x  1 0 1  <  Re*,  <  10^,  it  is  apparent  that  the  form  of  the  Cp  versus  Mro 
variation  at  0.  9  <  <  1.  1  is  a  function  of  free -stream  Reynolds  num¬ 

ber.  For  high  Reynolds  numbers  (i.  e. ,  Re*  >  105),  there  appears  to 
be  a  smooth  transition  from  subsonic  to  supersonic  drag  values.  For 
2  x  102  <  Rea  <  10^,  there  is  not  a  smooth  transition  from  subsonic  to 
supersonic  drag  values.  For  2  x  101  <  Re*,  <  2  x  102,  the  change 
from  subsonic  to  supersonic  drag  values  is  accompanied  by  relatively 
large  changes  in  CD  for  small  changes  in  M^. 

The  significance  of  the  9  results  is  that  they  indicate  a  basic  diffi¬ 
culty  in  measuring  sphere  drag  at  M*  =  1.  0  in  a  short  ballistic  range 
such  as  the  VKF  100-ft  Range  K.  The  reason  for  this  is  that  a  velocity 
drop  of  at  least  10  ft/sec  is  required  to  derive  an  accurate  "3rag  coeffi¬ 
cient  at  M,  =  1.0.  This  means  that  for  =  1.  0  the  Mach  number 

varies  from  0.  996  to  1. 004.  At  low  Reynolds  numbers  in  this  speed 
regime,  it  is  not  completely  valid  to  assume  that  Cp  is  constant,  and 
the  VKF  1000-ft  range  (Hyperballistic  Range  (G))  would  be  required  in 
order  to  accurately  detect  the  change  for  a  AV  of  only  10  ft /sec. 
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In  Fig,  7  the  slashed  line  in  the  region  of  M,,,  =  1,0  represents  the 
authors'  best  guess  as  to  the  value  of  Cjj.  In  applying  these  values  to 
the  full-scale  sphere,  a  reduced  degree  of  confidence  should  be  ascribed 
to  upper  air  density  values  obtained  when  Mffl  =  1.0. 

As  an  item  of  general  interest.  Appendix  III  presents  schlieren 
photographs  of  sphere  shock  wave  development  at  transonic  speeds. 

4.5  SUPERSONIC  SPHERE  DRAG  COEFFICIENTS 

A  review  of  Figs.  6k  through  v  reveals  that  there  is  a  limited  amount 
of  data  in  the  1.  0  <  <6.0  regime  from  other  sources  for  Ty^/T^  =  1.  0 

(Refs.  5,  6,  and  8  through  12).  These  data  provide  a  good  indication  of 
the  variation  of  drag  coefficient  with  Reynolds  number  at  a  fixed  Mach 
number  for  Re*,  >  104.  With  the  exception  of  some  data  obtained  at 
1.  45  <  M^,  <  1.  85  and  200  <  Re,,,  <  300,  these  data  are  also  in  good  agree¬ 
ment  with  present  data  for  Re,,  <  104. 

It  is  of  interest  to  consider  the  form  of  the  Cjy-Re,,,  curves  shown  in 
Figs.  6k  through  v.  The  total  drag  coefficient  of  a  sphere  can  be  written 
as  the  sum  of  drag  components  attributable  to  forebody  pressure,  after¬ 
body  pressure,  and  friction,  viz 

Cd  =  CDPF  4  CDPB  4  C[)p 

The  forebody  pressure  drag  (Copp)  of  a  sphere  can  be  derived  from 
the  pressure  distribution  over  a  hemisphere  (cf.  Refs.  23  and  24).  Clark 
(Ref.  24)  has  derived  a  simple  expression  for  forebody  drag 

CDPF  =  °-901  '  0.462/M^2 

This  derivation  is  in  good  agreement  with  Hoerner's  (Ref.  23).  From  a 
consideration  of  the  experiment;!!  measurements  of  stagnation  pressure 
on  source-shaped  bodies  by  Sht  nan  (Ref,  25)  shown  in  Fig.  11,  this 
pressure  term  would  be  expected  to  be  essentially  constant,  for 
Re^  >  3  x  102. 

Lehnert  (Ref.  26)  has  measured  the  afterbody  pressure  drag  of  a 
sphere  over  a  range  of  Reynolds  numbers.  Tnese  data  are  shown  in 
Fig.  12.  Some  afterbody  pressure  drag  values  have  been  derived  from 
these  curves  for  Re,,  -  109  and  are  reported  with  some  similar,  more 
recent  data  by  Jerrell  (Ref.  27)  in  Fig.  13.  The  two  sets  of  data  are  in 
good  agreement. 
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Also  shown  in  Fig.  13  is  the  summation  of  the  two  terms  Cj)pjr 
and  Cjp pg  compared  with  the  ballistic  range  values  of  total  drag.  The 
small  difference  between  these  two  curves  is  representative  of  the 
friction  drag  component  which  for  a  sphei'e  at  Re^  =  10®  would  be  ex¬ 
pected  to  be  small. 

At  supersonic  speeds  and  for  10^  <  Re,,  <  10®  the  sphere  drag 
coefficient  decreases  as  the  Reynolds  number  decreases  from  10®  to  lo\ 
e.  g. ,  Fig.  6o.  As  noted  earlier,  the  value  of  C£>pjr  is  substantially 
constant  over  this  range  of  Reynolds  numbers.  With  decreasing 
Reynolds  number  the  friction  drag  component,  Cp>p,  increases.  Thus, 
the  decrease  in  Cp-p  must  be  explained  in  terms  of  a  decreasing  value 
of  Copp.  This  is  in  agreement  with  Lehnert's  (Ref.  26,  Fig.  12) 
measurements  of  Cppp  in  the  wind  tunnel,  which  show  a  similar  varia¬ 
tion  with  Reynolds  number. 


4.6  EFFECT  OF  MODEL  WALL  TEMPERATURE  AND 

ACCELERATION  ON  SPHERE  DRAG 

In  determining  the  correct  drag  coefficient  to  use  in  the  supersonic - 
hypersonic  speed  regime,  attention  must  be  paid  to  the  wall  tempera¬ 
ture  of  the  sphere.  The  effect  of  wall  temperature  upon  the  viscous 
drag  of  a  sphere  at  near-continuum  flow  conditions  has  been  theo¬ 
retically  demonstrated  by  Davis  and  Flugge-Lotz  (Ref.  28).  Hayes  and 
Probstein  (Ref.  29)  present  equations  showing  the  effect  of  wall  tem¬ 
perature  on  drag  coefficient  for  free -molecule  flow  conditions.  On  the 
basis  of  this  information  and  several  experimental  studies,  it  is  obvious 
that  the  drag  coefficient  in  the  transition  regime  between  free -molecule 
and  continuum  flow  conditions  is  a  function  of  wall  temperature. 

Discussions  with  Messrs.  Wright  and  Morrissey  of  Cambridge 
Research  Laboratory,  who  are  interested  in  the  passive  and  instru¬ 
mented  spheres,  respectively,  have  indicated  that  they  expect  the  wall- 
to-free -stream  temperature  ratio  of  these  spheres  to  be  close  to  unity 
over  the  entire  trajectory.  Thus,  for  subscale  experimental  data  to 
cc. erectly  simulate  the  full-scale  event,  this  temperature  ratio  should 
be  near  unity.  With  few  exceptions,  none  of  the  available  supersonic 
wind  tunnel  data  obtained  to  date  has  fulfilled  this  requirement. 

/ 

The  ballistic  range  data  obtained  both  here  and  elsewhere  for 
T.  0  <  M„  <6.0  have  been  obtained  over  relatively  short  flight  distances, 
and  consequently  no  appreciable  model  heating  would  be  expected. 
Therefore,  it  has  been  assumed,  on  the  basis  of  approximate  calcula¬ 
tions  of  the  heating  rates,  that  the  ballistic  range  data  c-  "respond  to 
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the  condition  where  the  wall-to-free-stream  temperature  ratio  is  unity. 
When  the  ballistic  range  and  full-scale  spheres  operate  at  conditions 
where  Tw/Ta>  *1.0,  it  is  possible  to  use  the  ballistic  range  data 
directly  for  the  full-scale  drag  values.  To  determine  whether  small 
changes  in  wall  temperature  could  affect  the  drag  coefficient  signifi¬ 
cantly,  a  study  was  made  of  several  sets  of  data  (Refs.  5  and  30 
through  39),  where  Tw/T(i)  ^  1.0.  These  data  are  shown  in  Figs.  14a 
through  f.  The  data  are  plotted  with  Cp>  as  a  function  of  Re2  where 
Rei  is  considered  to  be  a  better  parameter  for  comparison  than  Re^ 
when  >  1.  Some  of  these  data  (cf.  Figs.  14e  and  f)  fall  outside  the 
Mach  number  limits  of  the  present  investigation  but  are  included  to 
illustrate  the  consistent  effect  of  wall  temperature  upon  sphere  drag 
coefficient. 

For  Figs.  14e  and  f,  the  Tw/Tu0  =1.0  data  have  been  derived  from 
VKF  ballistic  range  data  reported  in  Ref.  5,  together  with  some  hitherto 
unpublished  data  obtained  between  1905  and  1909.  Sphere  drag  data  for 
8.  5  <  M,,,  <  21.  5  are  snown  in  Figs.  15a  through  c,  At  low  Reynolds 
numbers  there  is  more  scatter  in  these  data  than  there  is  for'  the  present 
investigation.  These  data,  together  with  the  results  of  the  present  in¬ 
vestigation,  are  summarized  in  F’ig.  10.  It  is  from  Fig.  10  that  the 
^w/T^  =  1.0  data  shown  in  Figs.  14e  and  f  are  derived. 

Figure  14a  demonstrates  the  consistency  of  experimental  meas¬ 
urements  from  a  variety  of  sources.  The  exception  to  this  good  agree¬ 
ment  is  the  data  presented  in  tabular  form  in  Ref.  34  derived  from 
University  of  Minnesota  data  referenced  in  Kef.  1. 

All  of  the  data  (Figs.  14a  through  f)  show  that  at  a  fixed  Reynolds 
number  the  drag  coefficient  increases  as  Tw/ increases.  To  illustrate 
this  effect  more  clearly,  the  data  contained  in  Fig.  14  have  been  re¬ 
plotted  as  Cp)  versus  Tyj/T^  at  various  discrete  values  of  Rc2  in  Figs.  17a 
through  f.  The  variation  of  drag  coefficient  with  velocity  and  wall  tem¬ 
perature  has  been  defined  theoretically  (cf.  Ref.  29)  in  the  free -molecule 
limit,  and  the  results  are  summarized  in  Fig.  18.  As  the  free -molecule 
limit  is  approached,  the  form  of  the  variation  of  drag  coefficient  with 
wall  temperature  should  approach  that  which  exists  in  the  free-molecuie 
limit.  The  M,,,  =  2.  0  data  contained  in  Fig.  17a  is  the  most  complete 
set  of  data  from  which  to  establish  the  effect  of  wall  temperature.  It  is 
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evident  from  these  data  that  when  Cg)  >  0.  7  the  temperature 

effects  upon  drag  coefficient  approximate  those  which  are  predicted 
theoretically  in  the  free -molecule  limit.  Therefore,  the  curves  that 
have  been  faired  through  the  data  in  Figs.  17b  through  f  have  been 
based  on  the  assumption  that  when  Cp>  >0.7  CpjjrM  the  form  of  the  Cg> 
variation  with  Tw/Ta)  is  that  given  by  theory  in  the  free -molecule  limit. 

In  Figs.  14b,  c,  and  e  the  Cp>  versus  Reffi  variation  for  T^/T^  ~  1.0 
is  shown  as  a  solid  line  and  a  dashed  extrapolation.  The  dashed  posi¬ 
tion  of  these  curves  has  been  derived  from  data  obtained  for  TW/T0D  >  1.0 
shown  in  Figs.  17b,  c,  and  e. 

In  discussions  with  J.  Morrissey  of  AFCRL,  it  was  indicated  that 
there  was  an  interest  in  obtaining  some  estimates  of  sphere  drag  at 
Re^  ®  1.  0.  Using  the  data  contained  in  Figs.  14a  through  f  and  17a 
through  f,  some  engineering  extrapolations  outside  the  scope  of  the 
experimental  data  have  been  made  to  satisfy  this  requirement.  The 
result  is  given  in  Fig.  19.  Obviously,  there  is  an  increased  uncertainty 
in  these  values  of  drag  coefficient. 

Tne  results  that  have  been  discussed  in  this  section  have  been  ob¬ 
tained  in  (1)  ballistic  ranges  where  the  deceleration  may  be  Irgh,  e.  g.  , 
>100  g;  (2)  free  flight  in  wind  tunnels  where  the  acceleration  is  low, 

0  [  1  g] ;  and  (3)  models  mounted  on  a  balance  in  a  wind  tunnel  where 
there  is  no  acceleration.  The  effect  of  acceleration  on  the  drag  of 
spheres  and  cylinders  has  been  investigated  by  several  authors 
(Refs.  40  and  41.  )  Any  differences  in  the  measured  drag  coefficients 
discussed  here  can  be  shown  to  be  caused  by  differences  in  wall  tem¬ 
perature,  This  suggests  that  acceleration  (or  acceleration)  has  had  no 
effect  upon  the  sphere  drag  coefficient  for  the  levels  encountered  in  the 
ballistic  range  during  these  tests. 


4.7  EFFECT  OF  SURFACE  IRREGULARITIES  AND  MODEL  SCALE 
ON  DRAG  COEFFICIENT 

The  ROBIN  Sphere  has  six  dimples  in  its  surface  caused  by  the 
holddown  points  of  the  internally  mounted  radar  reflector.  An  approxi¬ 
mation  to  these  surface  irregularities  has  been  made  on  some  spheres 
launched  in  support  of  this  test  program.  Also,  an  investigation  has 
been  made  of  physical  model  scale  at  one  specific  condition.  The  results 
are  listed  below. 
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Those  results  indicate  that  for  the  conditions  of  this  investigation  there 
are  no  significant  effects  of  scale  or  surface  irregularity  on  sphere 
drag  coefficient  which  is  seen  to  vary  only  ±1  percent  from  the  mean 
in  all  cases  represented  in  the  above  table.  Sivier  (Ref.  17)  tested 
some  moderately  roughened  spheres  at  Re*  <  2  x  10^  and  found  no  effect 
upon  the  sphere  drag  coefficient  at  ~  0.  16. 

4.8  DRAG  COEFFICIENT  AT  HIGH  SPEEDS 

As  can  be  seen  from  Section  5.  6,  there  is  not  a  large  body  of  ex¬ 
perimental  data  concerning  the  effects  of  wall  temperature  upon  drag 
coefficient.  Consequently,  to  make  use  of  as  much  of  the  available 
data  as  possible  concerning  the  effect  of  wall  temperature,  it  has  been 
necessary  to  consider  data  at  speeds  higher  than  that  required  in  the 
present  investigation.  This  has  resulted  in  the  derivation  of  Cp>  versus 
Mm  curves  at  fixed  Reynolds  numbers  shown  in  Fig.  16. 

At  first  sight,  the  form  of  the  Cp>  versus  curve  at  low  Reynolds 
numbers  appears  to  exhibit  an  unexpected  variation.  However,  it  must 
be  remembered  that  Knudsen  number  is  a  more  realistic  indicator  of 
the  approach  io  the  free-molecule  limn  rather  than  Reynolds  number. 

In  fact,  Knudsen  number  is  proportional  to  M^/Re^.  This  means  that 
for'  Re^  =  20  the  Knudsen  number  at  =  20  is  twice  that  for  =  10. 
This  in  turn  means  that  the  =  20  condition  is  closer  to  the  free- 
molecule  limit  than  the  -  10  condition  even  when  Rea  is  the  same. 
This  explains  why  the  high  Mach  number  drag  values  are  closer  to  the 
free-molecule  limits  than  the  low  Mach  number',  values . 

SECTION  V 
CONCLUSIONS 

Measurements  have  been  made  such  that  reliable  values  of  sphere 
drag  coefficient  may  be  derived  for  any  value  of  and  Re^  for 
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Tw /Tg,  ~  1.0  within  the  bounds  0.  1  <  <6.0  and  2  x  10^  <  Re^,  <10^ 

with  an  uncertainty  no  larger  than  ±2  percent. 

Based  on  these  and  other  published  data,  there  is  sufficient  informa¬ 
tion  contained  herein  to  predict  the  effect  of  wall  temperature  on  Cp> 
when  4  1.0  for  2,  0  <  <6.0. 

From  the  present  results  and  those  obtained  elsewhere,  no  meas¬ 
urable  effect  of  acceleration  upon  drag  can  be  detected,  and,  therefore, 
no  effect  is  expected  under  the  conditions  of  the  ROBIN  and  "instrumented 
Spheres"  programs. 

A  brief  inveecigation  of  model  scale  and  surface  irregularities  indi¬ 
cates  that  for  Mro  =  0.  25  and  Re,,,  =  10^  no  effect  on  drag  could  be 
measured. 

There  is  reasonable  agreement  between  the  present  low  speed  data 
( <  0.  25  and  Re^  <  10^)  and  the  classical  data  which  has  resulted  in 
the  derivation  of  the  "standard  drag  curve.  "  Any  differences  in  the  two 
sets  of  data  may  be  explainable  in  terms  of  velocity  differences  in  the 
basic  data. 

The  present  investigation  represents  the  most  comprehensive  experi¬ 
mental  program  to  date  to  define  sphere  drag  in  the  velocity -altitude 
regime  of  interest  in  the  "falling  sphere"  programs  for  measuring  upper 
air  density.  It  is  believed  that  the  present  Cp>,  and  Re,,  map  is  the 

most  accurate  produced  to  this  time. 
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b.  20-mm-Cannon  Sabots 


d.  Typical  Spheres  Used  in  this  Test 

Fig.  3  3 
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a.  Cold-Gas  Gun  Sabots 


v  *  * 


c.  Rifle  Sabot 

Size  and  Material  of  the  Spheres,  from  left  to  right, 
respectively,  areas  follows: 


1-in.  -diam  Nylon 
1/2 -in.  -diam  Copper 
7/16-in.  -diam  Oylite 
3/8-in.  -diam  Stainless  Steel 
1/4-in.  -diam  Aluminum 
1/4-in.  -diam  Beryllium-Copper 
1/8 -in.  -diam  Lexan 
3/32 -in.  -diam  Nylon 


bots  and  Spheres 


29 


2 


3 


i 


1.  MODEL  AND  SABOT  TOGETHER  IN  LAUNCH  TUBE. 

2.  MODEL  AND  SABOT  SEPARATED  AFTER  SABOT  STOPPED  IN 
LAUNCH  TUBE. 

3.  SPHERE  GOES  ON  DOWN  RANGE. 


a.  Cold-Gas  Gun  Stripping 


12  3  4 


1.  MODEL  AND  SABOT  10GETHER  IN  LAUNCH  TUBE. 

2.  MODEL  AND  SABOT  SEPARATED  AFTER  PASSING  THROUGH 
PIN  OR  WASHER  STRIPPER. 

3.  SABOT  STRIKES  ANGLED  RAMP  AND  DEFLECTS  VERTICALLY. 

4.  SABOT  STRIKES  CATCHER  PLATE  -  SPHERE  PASSES  THROUGH 
HOLE  AND  ON  DOWN  RANGE. 

b.  20-mm-Cannon  Stripping 


1  2  3 

1.  MODEL  AND  SABOT  TOGETHER  IN  RIFLED  LAUNCH  TUBE. 

2.  PETALLED  SABOT  SPREADING  UNDER  ACTION  OF  CENTRIFUGAL 
FORCES. 

3.  SABOT  ARRESTED  BY  CATCHER  PLATE  -  SPHERE  PASSES 
THROUGH  HOLE  AND  ON  DOWN  RANGE. 


c.  Rifle  Stripping 

Fig.  4  Model  Separation  Techniques  Used  in  Range  K 
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SHOT  2627 


Z 


Fig.  5  Typical  Least-Squares  Linear  Fit  to  Velocity-Distance  Data 
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VKF,  Ballistic  Range,  1970 
Laurence,  Ballistic  Range,  Ref.  2 
Lawrence,  Ballistic  Range,  Unpublished 


°  VKF,  Ballistic  Range,  1970 
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Fig.  6  Continued 
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p.  1.65  <  <  1.85 


q.  1.9  <  IVL  <  2.2 
Fig.  6  Continued 
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b.  0.2  CM.  <1.6 
Fig.  7  Continued 
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Fig.  7  Continued 
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e.  1.0  <  M„  <  6.0 
Fig.  7  Concluded 
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Fig.  8  Comparison  of  Sphere  Drag  Measurements  at  Low  Speeds 
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Fig.  9  Variation  of  Sphere  Drag  Coefficient  with  Mach  Number  {Pre-1970  Data) 
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Fig.  11  Stagnation  Point  Prest  »re  as  a  Function  of  Mach  Number  and 
Reynolds  Number 


- Afterbody  Pressure  Drag,  Cqp0,  Ref.  26 

- Torebody  Pressure  Drag,  Cqpp,  Ref.  24 

CDPB  +  CDPF 

-  Total  Drag  -  Aeroballistic  Range,  cf  Fig.  7d 

o  Afterbody  Pressure  Drag,  Cqpb,  Ref.  27 

Rero  ~  106 


Fig.  13  Total  Drag  as  a  Function  of  (vlach  Number 
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V.KF.  Ballistic  Range 
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Fig.  14  Continued 
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Fig.  17  Variation  of  Sphere  Drag  Coefficient  with  Wall  Temperature 
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Fig.  19  Variation  of  Sphere  Drag  Coefficient  with  Mach  Number  -  An 
Extrapolation  to  Re„  =  2.0 
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TABLE  I 

SPHERE  MATERIALS  AND  SIZES  USED  IN  THE  RANGE  TEST 


Diameter,  in. 

3/32 

1/8 

1/4 

3/8 

7/16 

— 

1/2 

1 

Copper 

V 

V 

— - , , , ,  — 

V 

Beryllium -Copper 

V 

Steel 

V 

V 

V 

Aluminum 

V 

V 

Nylon 

V 

V 

si 

V 

V 

Lexan®* 

V 

Dylite®  ** 

V 

sj 

V 

*Lexan,  a  registered  trademark  of  the  General  Electric 
Company,  is  of  density  slightly  greater  'than  nylon. 

**Dyliie,  a  registered  trademark  of  Sinclair  Koppers  Company, 
is  an  expandable  styrofoam,. 
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TABLE  II  (Continued) 


APPENDIX  III 

DEVELOPMENT  OF  SHOCK  PATTERNS  AT  TRANSONIC  SPEEDS 


The  following  schlieren  photographs  are  included  in  this  report  as 
an  item  of  general  interest.  These  pictures  show  the  development  of 
bow  and  wake  shock  patterns  for  a  0.  25 -in.  -diam  nylon  sphere  traveling 
at  speeds  near  --  1.  For  these  conditions,  1050  <  Re,,,  <  1250. 


Flight  Direction 


t 

Bow  Shock 

Mrj)  * ).  012 


t 

Wake  Shock 


Fig.  Ill-*  Development  of  Shock  Patterns  at  Transonic  Speeds 
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Bow  Shock 

Wake  Shock 

Bow  Shock 

Wake  Shock 

o  *  1  045 

Mgj  •  1. 130 

Fig.  lli-1  Concluded 
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